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SUMMARY 

The anomalous Behavior of the thermal conductiv- 
ity of five gases, carhon dioxide, argon, nitrogen, 
oxygen, and methane has Been investigated using the 
theoretical approach of Brokaw and the equation of 
state given By Bender. The analysis presented herein 
properly descrihes the trends in the data and in gen- 
eral is in good agreement with the data for p/p(. < 
1.1. The theory falls 10 to 15 percent above the 
data at higher densities and appears about 10 percent 
high in the peripheral regions. The results seem ad- 
quate for engineering calculations and useful in es- 
tablishing the thermal conductivity Behavior of 
fluids for which no data have yet Been taken. 

NOMENCLATURE 

A,B constants in Eq. (15) 

Cp specific heat at constant pressure, 

j/gm-K 

Cp specific heat at "zero" pressure, j/ 

gm-K 

Cp =Cp-Cp^ "reacting" specific heat, j/gm-K 

Cp level adjustment constant 

Dp self diffusion coefficient, cm^/sec 

Dpp Binary diffusion coefficient, cm^/sec 

^ molecular weight, gm/gm-mole 

n number of monomers 

P pressure, MN/m^ 

E gas constant, j/gm-K 

T temperature, K 

xP parameter defined By Eq. (iz) 

Y parameter defined By Eq. (ll) 

Z compressibility 

A thermal conductivity, w/cm-sec 

Aj; contribution to thermal conductivity due 

to dissociation, W/cm-sec 



modified form of Ar, w/cm-seo 

Aj> 

frozen thermal conductivity, W/cm-sec 

a| 

thermal conductivity at "zero" density, 
w/cm-sec 


viscosity, gm/cm-sec 

f(p) 

a function of density, W/cm-sec 

P 

density, gm/cm^ 

(J 

diameter of the molecule. Angstroms 

sd.D* 

collision Integral - self diffusion 

a( 2 , 2 )* 

collision integral - viscosity 

>^i 

monomer in the cluster 


monomer 

Subscripts : 


c 

critical 

pseudo 

a representative value 

CO 2 

carbon dioxide 

% 

nitrogen 

AE 

argon 

CH 4 

methane 


oxygen 

INTRODUCTION 



During the past decade, the debate on the exist- 
ence of a "spike" in the thermal conductivity at the 
thermodynamic critical point has flourished. One of 
the earlier investigators to note this anomalous Be- 
havior was Guildner (l). Using a concentric cylinder 
cell, he found a marked increase in the thermal con- 
ductivity of carbon dioxide near its critical point. 
His conductivity data for selected Isotherms, when 
plotted versus density, are nearly symmetrical about 
the critical density with some additional Increase at 
the lower densities. Guildner (l) also suggested 
thermal conductivity to Be extremely large if not of 
infinite value at the critical point. 
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in Eq. (4) is 


Sengers (2, 5, 4) refined the experimental pro- 
cedure. Using a parallel plate cell he acquired some 
excellent thermal conductivity data in the near crit- 
ical region for carhon dioxide. He took great care 
to minimize cell convection and through careful anal- 
ysis was ahle to show it to he negllglhle. In gen- 
eral, Sengers' peak conductivities are less than 
Guildner's for equivalent Isotherms, perhaps indica- 
ting some undetected convection effects in Guildner's 
data. While some convection will prohahly he present 
in measurements .using cylindrical cells, Bailey (5) 
took near critical argon data in hoth the horizontal 
and vertical positions without any notieeahle gravi- 
tational effects . 


The "reacting" specific heat Cp^ 
defined as an excess specific heat 

Cpr = Cp(P,T) - Cpq(T) (5) 

and Gj_ is a level adjustment constant. This con- 
stant is determined hy matching the theory to a single 
experimental data point. Once C]_ is determined, the 
thermal conductivity values for the entire region can 
he found. 

Dp ‘is the self diffusion coefficient and is 
given as 


Although the existence of the "spike" has heen 
verified, procedures to predict this anomaly for 
various fluids are lacking. Herein we will combine 
the equation of state of Bender, Ref. 6, the theo- 
retical approach of Brokaw, Ref. 7, the excellent 
carhon dioxide data of Sengers, Ref. 2, 3, 4, and 
other data to formulate a technique to predict the 
anomalous behavior of the thermal conductivity in the 
near critical region, for reduced densities in the • 
range 0.6 < p/p^ < 1.5. 

It should he noted that Benders state equation 
Ref. 6, is analytic at the critical point. Thus in 
principle one should not expect to obtain a complete- 
ly rigorous mathematical discussions of the anomaly 
which appears to he a mathematical singularity (2). 
However this does not seem to he of great importance 
for engineering applications. 

THECRETICAL BACKSlOUin) 

Brokaw' s Analysis and Modifications 

Brokaw, Ref. 7, assumed that the near critical 
thermal conductivity could he described in terms of 
a reacting gas. While he recognized that there 
exists a spectrum of cluster sizes, for simplicity, 
the gas was assumed to consist of a monomer and a 
single large cluster of n - monomers, l.e. 

nXp Xn (l) 

Based on the success of Brokaw 's theoretical analysis 
in predicting the thermal conductivity of the 
N20472-2H02 system. Ref. 8, this technique seems to 
he a sound approach. While the theory is rather in- 
volved, Brokaw represents the thermal conductivity 
as. 


•D-, = 


* 2.6693xl0‘5*^^/(pa2ffi(l^l)) 


( 6 ) 


where 


pPj 

n 


6 
5 


n(i,D* 


■ « 1.32 


(7) 


can be used to find Dp. The ratio of the binary to 
self diffusion coefficient (Epn/Di) was argued hy 
Brokaw, using scaling rules, to he of the form 

T 2 y 


"In 


Ll + 


(2.l)nV^ 




( 8 ) 


where the constant 2.1 is related to the geometry and 
density of the cluster and determined hy Brokaw, 

Ref. 7, using SP/3pJip as determined hy Sengers to 
give a good fit to the experimental carbon dioxide 
data of Refs. 2, 3, 4. Herein, it was found most ex- 
pedient to retain the 2.1 value of Brokaw and adjust 
the level for each fluid through Cj_, (Eq. (4)). 

The. number of monomers in a cluster, n, was given 
hy Brokaw and slightly modified herein: 


RT 

.Sp 


/(I - Z) 


(9) 


The absolute value was used to expedite computing 
in the peripheral regions when the symmetry constraint 
was not used. More properly, n should approach 1 
and the reaction contribution approaches zero, in 
these regions. 


X = 7^ + Xp 


( 2 ) 


SIMMETRY CONSTRAUST 


where is the sm of the contribution due to dis- 

sociation and diffusion of the clusters, and Xj> is 
the rest; i.e. the thermal conductivity at "zero" 
density X^ and the density effect contribution 
where 

Xp = Xt + f(p) (3) 


We found hy observation, that the data of 
Sengers, Refs. 2, 3, 4, possessed a high degree of 
symmetry about p/p^ = 1 for reduced densities in 
the range 0.6 < p/pp <1.5. Guildner's earlier data 
(l) also reflected this symmetry. We, then, decided 
to neglect the slight asymmetry at lower density 
levels and for this prediction to define a pseudo 
value of density for 1. < p/p^ <1.5 as 


is a function of density and X^ depends only on 
temperature; the form of f(p) was taken from Ref. 9 
for the fluids considered. 

The reacting contribution as defined herein is 
Xr = Cx>R = CxpCpjJ)i(Din/Dx) (4) 


Ppseudo 


=^Pc 


( 10 ) 


And since the computations are along an Isotherm 
(isobar for argon) such a combination of (ppseudoy ^) 
gives a complete symmetry with respect to p/p^ = 1; 
thus Eq. (4) when evaluated at (ppgeudo^ yields 
symmetric results. For example, the values of Xj 
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a-t p/pf, = 1.2 with the symmetry constraint, 

Eq. (10), are Identical to the values of 7^ at 
p/p(, = 0.8 because 

^p£seudo^ = (0.8, T) 
SEFGERS-KEYES AHALYSIS 


Recently, Sengers and Keyes, Ref. 10, devised a 
scaling procedure for determining the near critical 
thermal conductivity of carbon dioxide. The essence 
of their work is illustrated in Fig. 1, taken from 
Ref. 10, where the data were correlated using, the 
parametric grouping: 


Y = 


(A - Ay) yp/pe 

■=; ® • 
3.05x10"= — 


I -.6 


( 11 ) 



Using Pig. 1 and Eqs. (ll) and (12), one finds 
that for XP <0.4, 


A - Ap 

3.05x10' 



f(p,T) 


(13) 


or thermal conductivity is a function of p alone. 
For xP > 3, 


3.05x10 


-.6 




(14) 


It was found for . 4 < xP ^3. the data of Fig. 1 
can be approximated (~3 percent error) by Eq. (15). 


A + 


(xP) 


B 

l/. 35 


(- 6 ) 


(15) 


where A » 1. and B » .9. Note that the exponents 
represent those used to parametrize Y and xP. 


ANALYSIS AND RESUITS 


behavior. Ref. 7, for selected isotherms. The 
305.23'f isotherm is compared to the data of Sengers, 
Refs. 2, 3, 4; as can be seen in Fig. 3, the data and 
theory agree quite well for p/p„ < .85 but the 
theory is about 18 percent high for p/pc l*^* 
important thing to note here is that while the pre- 
diction does not precisely match the data the general 
trends are in good agreement. 

In order to obtain better agreement between the 
data and the theoretical prediction, the level was 
adjusted' to fit the point at 305.23 K and p/p^ = 0.8. 
At the same time the symmetry condition, discussed in 
the Theoretical Section, was invoked and the constant 
in Eq. (4) was found to be 1.15, i.e.. 


In Fig. 3 the thermal conductivity for the 
305.23 K isotherm, computed using Eq. (16), is again 
compared to the data of Sengers, Refs. 2, 3, 4. The 
results again indicate good agreement for p/p„ < l.lj 
but the theory is about 10 percent high for p/pj, » 
1.4. 


In Fig. 4, the computed thermal conductivity 
using Eq. (16) can be compared to the data of Sengers, 
Refs. 2, 3, 4, for several selected isotherms. The 
data and theory are in good agreement for the 304.33 
isotherm but the theory is about 13 percent high for 
the 313.13 K isotherm near p/p;, « 1.4. While agree- 
ment seems good near the critical point, further de- 
parture seems to indicate that the theory will over- 
predict the data by 10 to 15 percent. However, in the 
absence of data or as a guide to prediction trends, or 
in engineering use, the error is tolerable. 

In Fig. 5. the thermal conductivity computed 
using Eqs. (ll), (12), and Fig. 7 are compared to the 
data of Sengers Refs. 2j 3, 4, for isotherms. The 
computed values, which are. in excellent agreement 
with that data as presented in Ref. 10, are herein 
15 to 20 percent above the data, see Fig. 5. For 
p/p = 1 and the 304.33 K isotherm, the difference is 
about 20 percent. One must note however, that Bender 
Ref. 6, lists p;, = 0.464 gm/cm^ and Tp = 304.21 K 
for the critical constants while Senger (2,3,4,10) 
lists p 2 = .467 gm/cm^ and T^, = 304.19 K. Using 
Sengers critical constants instead of Benders lowers 
the predicted value by 8 percent. This illustrates 
the sensitivity of the results to selection of the 
critical constants. 

Thus it appears that either technique, Keyes and 
Sengers (lO) or Brokaw (7) will adequately describe 
the near-critical thermal conductivity of carbon- 
dioxide. The advantage of the technique of Ref. 10 is 
that it does not require a good equation of state in 
order to use the theory. The disadvantage is its ex- 
treme sensitivity to the selection of critical con- 
stants. 


Carbon Dioxide 


Nitrogen 


The thermal conductivity of carbon dioxide has 
been carefully assessed by Sengers, Ref. 2, 3, 4. 
Fig. 2 illustrates the use of Bender's equation of 
state. Ref. 9, (for derived properties, see Ref. ll) 
and Brokaw' s theoretical prediction of the anomalous 


With the fit to the carbon dioxide data as per 
Eq. (16), the question arises as to how well the 
"theory" will work for other fluids. In Ref. 12, 
Sengers reanalyzed, the nitrogen data of Ziebland and 
Burton, Ref. 13, and based on two near critical data 


tDue to a 0.02 K difference in critical temperatures of Refs. 2 ahd 6 the magnitude of the carbon dioxide 
Isotherms of Ref. 2 were decreased by 0.02 K. This permits a more proper comparison between the data sources, - 
Ref. 2 for thermal conductivity and Refs. 6 and 11 for PVT and derived properties. 
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points estimated the thermal conductivity curve for 
the 133,15 K isotherm. Using the theory of Brokaw, 
Ref, 1 , as modified herein and Bender's equation of 
state, Ref. 6, the level was adjusted to fit the 
point as p/p,, = 0,97 and T = 133.15 K. The con- 
stant of Eq. (4) was found to he 1.35, i.e., 

Alj-g = Tvp + 1.35 Ar (17) 

The computation was then completed using Eq. (17) 
along the entire 133.15 isotherm and was found to 
match the other data point at p/p„ » 1.49. See 
Fig. 6. , The estimated thermal conductivity curve 
given hy Sengers, Ref. 12, and the computed curve are 
in good agreement except for p/p^, *>• 0. 6 where the 
theory is about 15 percent high. Several other iso- 
therms were computed using Eq. (4), both below and 
above the critical pressure. These results are illu- 
strated in Pig. 6. 

The unusual behavior of the predicted thermal 
conductivity beyond 1.5 p/p^, is believed to be due 
to errors in the derivative SP/Sp and an overex- 
tension of the theory. 

Argon 

In Refs. 10, 11, Bailey presented data for the 
thermal conductivity of Argon for several isobars, 

(as opposed to isotherms). The level was again ad- 
justed to fit the point at p/p^ « .65 at 49.8 kg/ 
cm^ (4.8837 MR/m^) to determine the constant of 
Eq. (4). Cl was found to be 1.2, i.e., 

^AR = Ap + 1. 2 (18) 

As can be seen in Pig. 7, there is a distinct shift 
in the data of Bailey, Refs. 5, 14, away from p/pj, = 
1 to a peak near p/p^ » 1.16. This shift is unex- 
plained. However there also does appear to be a 
small shift in the peak Am as computed from 
Eq. (4) for the 50.8 kg/cm^ (4.9817 Mll/m^), see 
Pig, 7. Since the level predicted by the theory and 
the data are in good agreement in other regions, 
further modifications were deemed unwarranted in view 
of the unexplained shift in the thermal conductivity 
data toward higher p/pj, values. 

The values of thermal conductivity were then 
calculated for several selected isotherms as shown 
in Fig. 8, 

Methane 

The methane thermal conductivity data of 
Sokolova and Golubev, Ref. 15, was most difficult to 
analyze. .Moreover, the methane equation of state as 
given by Bender, Ref. 6, is probably the least well 
behaved, in the near-critical region, of the five 
gases correlated by Bender. In this case, the level 
was not adjusted to fit any particular point, but to 
best represent the trends of the 210 K and 216.3 K 
Isotherms. As can be seen in Pig. 9, the 210 K iso- 
therm data peaks near p/pj, » 0.7. There seems to be 
no apparent reason for this peak to occur as it does 
and is, as of now, unexplained. The 216.3 K iso- 
therm data possess the more conventional appearance 
and are in better agreement with the trends of both 
the theory and other data, e.g. see Pig. 4. The con- 
stant C]_ of Eq. (4) was selected as 2.2. 

^H4 = Ap + 2.2 Ar (19) 


It is apparent from Fig. 9, that either the data 
of Ref. 15 for the 194.85 K isotherm are too high for 
the p/p(j 1 region or the derived properties. 

Ref. 11, are too small. Also the adjusting constant 
(2.2) of Eq. (4) is much greater than for the other 
gases,' As such it also seems out of line. < 

Assuming that Eq. (19) does represent the best 
fit to these data, the conductivity of methane can be 
computed. See Fig. 10. In Pig. 10, the dip in the 
250 K isotherm can best be explained by the signifi- 
cant departure from the near critical region, i.e. an 
overextension of the theory, and troubles in the de- 
rived properties. Ref. 11. The Ap or frozen thermal 
conductivity for the 250 K isotherm is included for 
comparison; it would seem that one is sufficiently 
far from the critical point such that Eq. (19) should 
no longer be applicable. 

Oxygen 

The equation of state for oxygen. Ref. 6, is 
quite well behaved and agrees very well with the re- 
sults of Weber, Ref. 16, see also Ref. 11. As a re- 
sult the theory should produce a reasonably good pre- 
diction of the thermal conductivity of oxygen. Assum- 
ing the level adjustment to be the same as for nitro- 
gen, 

^Og = ^p + 7^ (20) 

the thermal conductivity of oxygen was computed for 
selected isotherms as illustrated in Pig. 11. 

LEVEL ADJUSTMENT CONSTAMT 

While an a priori technique to predict the value 
of C^, Eq. (4), has not yet been established, there 
does appear to be a trend with molecular weight as 
seen in Pig, 12. This tendency to deviate with mol- 
ecular weight may be related to geometry and density 
of the monomer which is in turn related to the con- 
stant in Eq. (8) and the number of monomer units (n), 
Eq. (9). If this trend is correct, a value of C]_ = 
1.3 should be used for oxygen. In the absence of 
data, 1.35 was used. Note that methane does not 
appear to follow the trend established by the other 
gases. 

CONCLUSIONS 

The theoretical approach of Brokaw combined with 
an adequate equation of state as described by Bender 
can be used to predict the anomalous behavior of the 
thermal conductivity in the near critical region for 
five fluids: carbon dioxide, argon, nitrogen, oxygen, 

and methane. In general the theoretical prediction as 
described herein is 10 to 15 percent higher than the 
data for p/p^, » 1.4. The theory for p/p,, <1.1 is 
generally in good agreement with the data. Methane 
is by far the most questionable of the theoretical 
predictions. 

In this analysis the level adjustment constant, 
Cj_, is determined by matching the theory at a single 
data point. Once Cj_ is determined, theimal con- 
ductivity values for the entire region can be found. 

A trend has been noticed for the level adjusting 
constant with molecular weight. This effect is prob- 
ably coupled with Brokaw' s empirical determination of 
a constant which is related to the geometry, density. 
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and tlie number of monomer units. 

The Keyes and Sengers technique also predicts . 
the anomaly. The advantage of their technique is 
that it does not require an accurate equation of 
statej the disadvantage is that it is extremely sen- 
sitive to the choice. of the critical constants. 

It is felt that these approaches can be used to 
predict trends and give satisfactory engineering re- 
sults. 
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Figure 1. - Correlation of the thermal conductivity of 
carbon dioxide (taken from ref. 6). 




p/Pc 

Figure 3. - Computed thermal conductivity of 
carbon dioxide adjusted for levei at p/p^ = 
0.8 and symmetry in for the 305. 23 K 
isotherm. Data from references 2 to 4. 


ISOTHERMS, 


o 304.33 
□ 305.23 



T = 304.21 K 
Pj-= 7.3835MN/m^ 
p- = 0. 464 GM/CC 

c I I 


.21 ^ ^ 1 1 

.6 .8 1.0 1.2 1.4 

P/pc 

Figure 4. - Comparison of computed thermal 
conductivity of carbon dioxide with the data 
of Sengers (refs. 2 to 4) for seiected iso- 
therms. The "frozen" thermal conductivity 
is included for comparison (ref. 9). 
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o 49. 8 KG/CIVl^ (4. 8837 AAN/m^) 
o VERTICAL Cai 'l50.8KG/CM^ 

A HORIZONTAL CELL J (4. 98177 MN/m^) 


WATT 
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Figure 7. - Comparison of computed thermal conductivity 
of argon with the data of Baiiey (refs. 5 and 14) for two 
isobars (P^. = 48. 392 Bailey, Pj. = 48. 014 atm PROGRAM). 
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Figure 8. - Computed thermal conductivity for 
argon for a few selected isotherms. 
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Figure 9. - Comparison of computed thermal 
conductivity of methane with the data of 
reference 15 for three isotherms. 
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Figure 10. - Computed thermal conductivity 
for methane for selected isotherms. 
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Figure 11. - Computed thermal conductivity of 
oxygen for seiected isotherms. 



MOLECULAR WEIGHT 


Figure 12. -Trend in the ievel adjustment 
constant Cj with molecular weight. 
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